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I. INTRODUCTION 



Giant resonances can be described as coherent superpositions of particle-hole excitations 
and their particle decays into given hole channels provide insight about the correspond- 
ing particle-hole amplitudes. Thus, the study of decay properties is a unique way to test 
model predictions for the resonance wave function. Experimental studies have recently been 
made of the particle decay properties of giant resonances by means of exclusive measure- 
ments where inelastically scattered projectiles, or light products of transfer reactions, are 
detected in coincidence with nucleons emitted by giant resonances of the target system. The 
corresponding results provide severe tests for the theoretical models. 

The Gamow- Teller resonance (GTR) is the manifestation of the spin-isospin nuclear 
zero sound. This mode was predicted theoretically as the mechanism for explaining the 
missing strength in the decay of /9-radioactive nuclei Its subsequent observation in (p,n) 
reactions has opened an important chapter of the historical development of experimental 
study of giant resonances when the availability of higher energy projectiles in the early 80's 
led to the realization of reaction conditions under which spin-isospin modes can optimally be 
excited 0. A first experiment aimed to study the proton decay of the GTR was performed 
in Groningen, and reported in Ref . || . The results of this measurement raised a number of 
questions, as they seemed to imply a vanishing spreading width for the GTR, which seems 
quite unlikely. 

Renewed interest in the study of the properties of the GTR is indeed testified by two 
recent measurements, with different characteristics, at Osaka and MSU. The proton decay 
of the resonance in 208 Bi has been investigated, and in the Osaka experiment the results 
appear in modest agreement with previous theoretical estimates given in Refs. || and ||, 
whereas preliminary indications from the MSU experiment seem to give a somewhat different 
result 0. This has motivated the present work, in which we apply to the charge-exchange 
modes a model of giant resonances which includes in a self-consistent way the coupling 
mechanisms leading to the damping of these modes. This model has been already used to 
study the properties of the isoscalar giant monopole resonance in 208 Pb ||. 

We have first concentrated, as a test case, on the isobaric analog resonance (IAR) in 
the nucleus 208 Bi where the data provides a clear picture of the resonance properties. The 
IAR would be degenerate with its isobaric multiplet partner, the ground state of 208 Pb, if 
the nuclear Hamiltonian were completely charge- independent. The energy difference, due 
mainly to Coulomb effects, has been studied over a long period of time and it is accurately 
known. The width of the IAR is known to be very small. This is because it has the same 
isospin as the parent state while surrounding states have the isospin of the ground state of 
208 Bi, which differs by one unit. Consequently, they couple only weakly with the IAR. Thus, 
excitation energy and width of this resonance provide a stringent test for any theoretical 
model of charge-exchange collective states. 

We shall show how this test has been passed by our model in Sec. |V|, while the predictions 
associated with the GTR are presented in Sec. |V|. Conclusions are drawn in Sec. [VT]. Before 
this, we discuss in Sec. [TI]the general features of the model already applied to other collective 
modes f| while more specific considerations relevant to charge-exchange modes are made 
in Sec. |ITI| . Detailed expressions for the Hamiltonian matrix and decay branching ratios are 
given in Appendices A and B. 
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II. GENERAL FORMALISM 



The Random Phase Approximation (RPA) provides, as a rule, an accurate description of 
the centroid of giant resonances and the fraction of the energy weighted sum rule (EWSR) 
exhausted by the mode. However, it is necessary to go beyond this approximation scheme in 
order to explain the damping properties of the collective motion. Indeed, giant resonances 
are known experimentally to have an energy width and therefore, a finite lifetime. A general 
theory of the resonance width can be found, e.g., in Refs. [j9 , |10|1 , and a list of references of 



much of the theoretical work carried out during the last decade can be found in ||11|| . In 
many of these calculations parameters were adjusted at various steps of the computational 
procedure in order to fit experimental data. In what follows we present a model in which a 
phenomenological effective interaction of the Skyrme type is used as the starting point but 
after this no further tuning of the parameters is introduced up to final results. 

We start by solving the Hartree-Fock (HF) set of equations for a given nucleus (^4, Z) 
using a Skyrme two-body interaction |fT2| . A set of occupied single-particle levels is then 
obtained, as well as a self-consistent mean field, which is then diagonalized on a basis made 
up with harmonic oscillator wave functions. This diagonalization procedure can also be 
replaced by solving the HF mean field with a box boundary condition. It provides a discrete 
set of occupied and unoccupied levels, many of which are at positive energy. The quasi- 
bound levels belong to this set. After choosing a convenient cutoff - details about numerical 



procedures are given in Sec. |T| - a finite set of occupied and unoccupied levels labeled by \i) 
is determined. In the following, we denote the corresponding energies and wave functions 
by Si and <pi (the radial part of the wave function will be expressed as Ui{r)). 

Let us call Qi the subspace of nuclear configurations made up with the HF ground state 
and all the possible one particle-one hole (lp-lh) excitations built within the set We 
shall denote by the same symbol subspaces and projectors onto them. Now, the nuclear 
Hamiltonian can be written as 

H = H + V ph , (2.1) 

where Hq is the HF Hamiltonian and V p h is the particle-hole (ph) interaction determined as 
the functional derivative of the self-consistent mean field with respect to the density. The 
usual Tamm-Dancoff approximation (TDA) or RPA in a discrete particle-hole space amount 
to solve TDA or RPA equations using the Hamiltonian QiHQi. 

As said above, giant resonances are known to have a decay width, and we can distinguish 
two main mechanisms which give rise to it. The energy of the vibrational nuclear motion 
can be transferred out of the system by an escaping nucleonB, or can be distributed among 
internal degrees of freedom giving rise to more complicated configurations than the initial 
configuration. The contributions to the width coming from these two damping effects are 
usually called respectively escape (r^) and spreading (r^) width. We must also mention 
another source of broadening of the line width, namely the Landau spreading which already 
appears at the level of discrete RPA or TDA calculations. 



Also a 7-ray can escape when selection rules allow it to do. In what follows we neglect this decay 
channel as it gives a minor contribution to the total width of the mode. 
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In order to account for escape and spreading effects, we build two other subspaces P 
and Q 2 - The space P is made up with particle- hole configurations where the particle is in 
an unbound state, orthogonal to all states To determine these unbound states, we use 
the following procedure. At positive energy e we solve, for each partial wave c = (I, j), the 
radial scattering equation for Hq, projected on the orthogonal complement of set We 
can thus ensure that the resulting outgoing wave functions uy~J have no overlap with any 
of the states \i), i.e., P and Q\ are orthogonal subspaces (for details, see [^,|TJ,0). The 
states have no resonant behaviour since quasi-bound states are among the states 

The space Q 2 is built with a set of "doorway states" , the first step in the coupling of the 
ordered resonance motion with the compound nuclear states, in which energy is distributed 
among all degrees of freedom in a statistical way. We denote these "doorway states" by 
\N). As a first approximation we can think of them as formed by 2p-2h configurations, but 
we discuss later a more physical choice of them in terms of states made up with one of the 
lp-lh excitations coupled to a collective vibration. 

We decompose the nuclear Green's function G as a sum of terms like QiGQi+QiGP+- ■ •. 
Using the technique described in |T(J one can show that by virtue of the equations of motion 
in the different subspaces and the properties of projectors, the Green's function QiGQi obeys 
an equation where the effective Hamiltonian, after truncation of higher order couplings, is 

U (u) = Q 1 HQ 1 + W\uo) + W l {uo) 

= Q 1 HQ 1 + Q X HP 1 , PHQ 1 + Q 1 HQ 2 \ , Q 2 HQ 1 , 

uj — PHP + ie uj — Q2HQ2 + ie 

(2.2) 

where u is the excitation energy. This energy- dependent, complex Hamitonian allows to 
work inside the space Q\. It has complex eigenvalues whose imaginary parts originate from 
coupling to unbound and to more complicated configurations and give rise to escape and 
spreading widths. 

The escape term W\uj) can be more easily evaluated if one replaces the complete Hamil- 
tonian H by the one-body part Hq. The neglect of matrix elements of Q\V p hP should be, 
in the present case, rather safe since discrete and continuum wave functions are essentially 
restricted to different radial intervals while V p h interaction has zero range. In this approxi- 
mation, the escape term can be written as 

W\uj) = lo-H -K, (2.3) 

where K is the inverse inside subspace Q\ of the Green's function containing only the 
mean field Hamiltonian. This Green's function can be easily computed and inverted as it 
has only one-body matrix elements whereas taking into account the two-body interaction 
would have resulted in a much harder task. Detailed expressions are given in Appendix 
A. The accuracy of the procedure has been checked in where the results obtained for 
the isoscalar monopole strength distribution in 40 Ca, by using the approximate of ([2. 3D, 
were compared with those of exact continuum-RPA calculations and excellent agreement 



between the two calculations was found (see Fig. 1 of [16|). Yoshida and Adachi made a 



similar comparison in [ffl and they concluded also that a good agreement can be obtained, 



provided a sufficiently large basis is employed (see Figs. 1-3 of [|TJ]). 
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The matrix elements on a basis of Q\ can be also determined in a straightforward way 
for the spreading term W^(u). We make the ansatz that the configurations \N) of Q2 are 
not interacting)!. This is reasonable since coupling among doorway states will correct the 
coupling between lp-lh states and doorway states in higher order of the residual interaction. 
If we take a basis of ph configurations, i.e., a set \ph) (particle levels are labeled by p,p' . . ., 
and hole levels by h, h! . . .), then the generic matrix element turns out to be 

N LU — U>N 

where we have omitted the subscript ph of the interaction, and uj^ is the energy of the state 
\N), in mean field approximation. Further clarifications about the evaluation of (|2.4|) are 



given in Sec. [II I 




It is found convenient, in order to solve the effective Hamiltonian ( |2.2| ), to work on the 
basis of the RPA or TDA states obtained by diagonalizing its first term. Using this basis, 
we provide explicit formulas for the matrix elements of the escape and spreading terms in 
Appendix A. In this matrix form, the eigenvalue equation for the effective Hamiltonian ( |2.2| ) 
is 

where T> is a diagonal matrix with the RPA eigenvalues, and the Ai matrices which are 
energy dependent contain the escape and spreading contributions, as explained in Appendix 
A. We have omitted the energy dependence of eigenvalues and eigenvectors in ( |2.5| ), for sake 
of simplicity in the notation, and we have collectively denoted respectively with F and F 
the amplitudes corresponding to positive and negative RPA eigenstates. The solutions of 
( |2.5|) are denoted by whereas the RPA basis vectors are labelled by \n), according to 

W = £^V>- (2-6) 

n 

It is worthwhile at this point to notice that an approach based on TDA instead of RPA 
leads to an equation in which only the upper-left quarter of the matrix that appears in ( |2.5| ) 
is present. The matrix which appears in (|2.5|) is complex symmetric as can be seen from its 
explicit form in Appendix A (see also ||13|| ). The transformation which makes it diagonal, 



that is the matrix of its eigenvectors, is complex orthogonal, that is 

F T F = FF T = 1. (2.7) 

A useful quantity characterizing a mode excited by a given operator O is the response 
function, 



2 These configurations display a rather high level-density in medium-heavy nuclei, of the order of 
10 2 -10 3 levels per MeV in the region of giant resonances if we think in terms of 2p-2h configurations 
[17], so that to calculate their mutual interaction would be a formidable task. 
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R(u) = (0|O + - 



H(u) + it 



The corresponding strength function is related to 



O|0>. 

by the well known relation 



1 



7T 



ImR(u) 



In terms of solutions of (|2.5| ) the strength function is 



S(u>) 



TV 



Im^2(0\O\u)'' 



1 



(2.8) 



(2.9) 



(2.10) 



where the squared matrix element of O appears, instead of its squared modulus, due to 
the properties of the eigenvectors which form a biorthogonal basis. We stress at this 
point that, although not explicitly written, the eigenvalues labeled by v as well as the 
wave functions needed to compute the matrix elements of O depend on u, as the effective 
Hamiltonian does. So, better than the eigenvalue distribution, which can change when 
diagonalizing the effective Hamiltonian at different energies, the strength function really 
carries information about the couplings taken into account, and can be directly compared 
with experimental results. Note also that, when all the above scheme is carried out without 
Q2 space, one should recover continuum RPA results. It is this kind of consistency check 
that we alluded to in mentioning Figs. 1-3 of and Fig. 1 of 



Another quantity which can be extracted from the model and which is actually measured 
in the particle decay experiments is the branching ratio B c corresponding to a particular 
decay channel. An escaping nucleon with energy e leaves a residual (A — 1) system in a hole 
state such that, by energy conservation Eh = e — u, where u is the initial excitation energy. 
The cross section a c for this decay was calculated in |T3j (see Eq. (3.40)), and a discussion is 



also given in [18] where a plane wave Born approximation (PWBA) is used to describe the 
reaction mechanism. In Appendix B we summarize this procedure and show how the same 
PWBA can lead from Eq. (3.38) of |13j to a more computable expression for the excitation 



cross section a exc . Using the results of Appendix B, the branching ratio comes out as 



B r . 



7) 



+ f(rv + r t 



where S uu > is given in (|B6| ) and 



dQ k j UtC (k)j*, c (k), 



with 



(2.11) 



(2.12) 



(2.13) 



In Eq. ( [2.13] ), ip c is the wave function describing the residual (A — I) nucleus in channel c, and 
(k) is the escaping particle wave function belonging to P space. It is shown in Appendix 



u 



(-) 



B that B c of Eq. ( 2.11 ) is a real quantity. The strength functions and the branching ratios 
are the main quantities we are going to show as results of the present model and which 
can be directly confronted with experiment. Before coming to this, we add some important 
remarks relevant to calculations of charge-exchange resonances. 
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III. CALCULATION OF CHARGE-EXCHANGE RESONANCES 



The model we have just presented was used in its TDA version and without the inclusion 
of the spreading term to calculate the IAR and GTR in 208 Bi || . Details about the HF 
and TDA procedures as well as the coupling with continuum configurations can be found 
in ||. On the other hand, the inclusion of the spreading term is an improvement of the 
model. Moreover, this term was never included in previous calculations of the IAR in the 
form of coupling with collective vibrations, as it is done here. While coupling to collective 
vibrations has been used for the GTR (cf. [^9j , p0[| ), the IAR case has only been treated by 
using for \N) the 2p-2h configurations to build ( JTH , [T7 1 ) . In the present work, we want 



to use the same model to construct the term for IAR and GTR, and therefore isospin 
conservation rules must be carefully treated as discussed at the end of this Sec. 

We have consistently used the two interactions SIII and SGII throughout the whole 
procedure. The HF mean field for 208 Pb is diagonalized on a basis of 15 shells of the set of 
harmonic oscillator eigenfunctions with Hujq=6.2 MeV (see also [0]). For the construction 
of subspace Qi, the set \i) contains all occupied levels and An unoccupied proton levels for 
each value of (/, j). An is 6 in the IAR case, and 7 in the GTR case, as in ||. The difference 
with this previous work is the inclusion in the present case of the Coulomb exchange term 
in the Slater approximation inside the HF mean field. 

We have then performed a TDA calculation within Qi, as the large neutron excess 
of 208 Pb should hinder ground-state correlations effects for proton particle-neutron hole 
configurations. This gives the basis of states |n) on which Eq. ( |2.5|) is solved. We have 
built the configurations of subspace Q2 by coupling a proton particle and a neutron hole in 
the discrete states of the set \i) with collective vibrations of the 208 Pb core. These doorway 
states, introduced above and labeled by \N) play an essential role in the damping of giant 
resonances. As pointed out in || , the reason is that the choice of 2p-2h as doorway states 
amounts to neglect a number of correlations induced by the residual interaction V p h- Many 
of these correlations are included if the low-lying collective vibrations are explicitly taken 
into account as doorway states. 

The collective vibrations of the nucleus 208 Pb have been calculated in self-consistent 
RPA, with the same SIII or SGII interaction used throughout the whole work, in a way 



completely analogous to P2 , where the coupling with collective vibrations was used to study 



the energies and spectroscopic factors of single-particle states. We refer then for details to 



22[, summarizing here the main points. We have built the complete spectra of isoscalar 
modes of J n equal to 2 + , 3~ and 4 + . In the configurations of Q2 however, only states with 
more than 1% of the total strength, and with energy less than 20 MeV, were included. 
The low-lying 5~ state was also included, but its contribution was found negligible. The 



collective states included in the calculations with SIII have been shown in 22 , and those 



corresponding to SGII are shown in Table [I] and compared with [|23|| . The limits in the 
accuracy with which this Skyrme-RPA calculation can reproduce the experimental findings 
are comparable to the case of similar calculations performed in the past decades. At this 
point we can build the matrix elements of the operator 

WV)=E V]W (3 . 1} 

N LU — U>N 



7 



as in Eq. ( |2.4j) (these are essential ingredients to build the matrix elements on TDA basis). 
^php'h' °f (H-ID is a sum °f terms whose diagrammatic representation is shown in Fig. |]. 
To evaluate these diagrams an expression for the interaction vertices is needed. This is 
provided by the particle-vibration coupling model, by introducing a one-body field which 
can be written in the form 

V = Y. E (a\tt\r)v(r)Y LM (r)\P) *W (3- 2 ) 

a/3 LnM 

In this equation, we have introduced the radial transition density Q^'{r) of the \n) states of 
the spectrum of phonons with angular momentum L. This is defined in terms of the RPA 
states as 

(n|^(r)^(r)|0) = Q^{r)Y LM {r), (3.3) 

where ip>(f) and if)(r) are the creation and annihilation operators of a nucleon at point r, 
respectively. The form factor v(r) appearing in ( |3.2| ) is related to the ph interaction derived 



from the Skyrme force by V^/ l (r 1 ,r 2 ) = f(r 1 )5(r 1 — r 2 ), as in p2 |. Once all these starting 
points are given, the evaluation of the diagrams of Fig. [I] is straightforward and their detailed 
expressions are given in Appendix A. 

Up to here, the procedure is the one already applied for the giant monopole resonance 
in ||. Since we deal here with charge-exchange modes, it must be noted that the operator 
can mix states with different isospins. The coupling ( |3.2|) is manifestly a scalar in the 
total fermion-boson isospin space. But the intermediate states \N) do not have pure isospin, 
as they contain a proton particle and a neutron hole (see the diagrams of Fig. |], where the 
phonons are isoscalar and therefore do not cause any change in the fermion isospin in the 
intermediate states, with respect to the initial state). This leads naturally to a coupling of 
the IAR, which has isospin quantum numbers |T, T z ) = |T ,T ) where T = — ^ ; with 
states which in general have a mixture of different T components. The nuclear part of the 
Hamiltonian actually forbids this coupling, and we impose that it is strictly forbidden (this 
amounts to neglecting Coulomb effects in the residual interaction). This can be done by 
projecting out the To component of the intermediate state \N). More precisely, we can write 
in isospin space 

\N) = c- 1 \N; T - 1, T -1} + Co| AT; T , T - 1) + c +1 \N; T + 1,T -1), (3.4) 

and we define three projectors P_i, P and P + i such that 

P i \N) = c i \N]T + i,To), — 1 < i < 1. (3.5) 

The coefficients q are simply geometrical factors (Clebsch-Gordan coefficients), 

c_i= +(2To-l)5(2T + l)-5, 

c = -(T + l)-5, 
c+x = +(T + 1)-5(2Tq + 1)-3. (3.6) 

Therefore, in the IAR case the correct spreading operator having the right isospin structure 
is 



S 



E 

JV 



VP \N){N\P V 



(3.7) 



Actually, the matrix elements of this operator are simply proportional to those of (|3.1| ), the 
multiplicative constant being |co| 2 = (T + l) _1 . In the GTR case, we use the operator x 
defined in a fully analogous way, and thus proportional to (|3.1|) by means of the constant 
|c~i| 2 . 

This procedure allows us to respect isospin symmetry. A similar procedure was not 



necessary in the IAR calculations of |15] and |17], where the intermediate states \N) were 
chosen as 2p-2h configurations. The interaction V was the nuclear two-body interaction, 
and the approximation used corresponds to the well-known Second RPA (SRPA). This 
automatically preserves the isospin structure, due to the respect of conservation laws in 
p4| , and to the character of the interaction. On the other hand, an approach similar 

for the GTR without any isospin factor, i.e., with the 
_ 1; differs by a factor |c_i| 2 close to 1 in the 



SRPA 



20 



to ours has been used in (TJ 
spreading operator ( |3.1|) ; as the correct one, 

case of 208 Pb, the numerical results are not markedly affected by this neglecting of isospin 
symmetry. 



w^- 



iv. ISOBARIC ANALOG RESONANCE RESULTS 

The calculated strength distributions in the IAR region are shown in Figs. (interaction 
SIII) and|3| (interaction SGII). Table |I| includes the main quantities which can be extracted 
from these strength distributions, that is centroid energy, width and percentage of strength 
of the resonance. From the TDA calculation, a discrete state emerges which, in SIII and 
SGII case respectively, lies at 18.59 MeV or 18.54 MeV, and carries 86% or 87% of the total 
strength (N — Z) of the operator T_ = Y,f=i t~(i). In fact, isospin is explicitly projected out 
only in the spreading term, whereas the TDA calculation is performed in a space with all 
isospin components included. Nevertheless, the strength relative to T_ operator selects as 
most collective state a single one, which was checked to be almost pure in isospin by looking 
at its T 2 expectation value. In Figs. and |3] the dashed line corresponds to a calculation with 
W* only, and the full line to the complete calculation. We remark first that, in contrast with 
other modes where a certain amount of fragmentation is present, in this case this strength 
function is essentially due to a single state. The effect of + is to shift slightly the 
peak energy and to produce a total width. With SIII, the peak is at 18.49 MeV and has a 
width of 152 keV whereas the corresponding values are 18.64 MeV and 99 keV with SGII. 
In both cases the state exhausts 97% of T_ strength. If one performs calculations with 
or only, one can obtain separately and TK We find that the sum + is equal, 
within a few percent, to the total width in the complete calculation. It can be seen that 
this width is rather sensitive to the interaction used, since the width calculated with SIII is 
50% larger than that of SGII. 

The complex collective state, expanded on the TDA basis, shows a squared overlap of 
0.90 (resp. 0.91) with the collective TDA state if one calculates with SIII (resp. SGII). 
We must also mention that the spreading term has been calculated with an averaging 
parameter iA added to the denominator of ( |3.7| ), for reason of convenience. This means 
that one makes a Lorentzian averaging of the distribution of intermediate states \N), with 
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A corresponding to half of the Lorentzian width. The quantity A has been chosen to be 100 
keV, but results are stable with respect to its variation, as doubling A implies an increase 
by only about 20 keV in the total width. Experimental values for IAR excitation energy 
and total width are respectively 18.8 MeV and 232 keV [p5] . This states exhausts more or 
less the whole T_ strength. 

Branching ratios corresponding to proton emission leaving the residual nucleus in a 
valence neutron hole state of 207 Pb are shown in Tables |TJ and [TV]. In these Tables, the 
column labelled "only W^" refers to a calculation with only the escape included, and the 
results correspond to those of 0. The other three columns under U W' + W^" include the 
results of the complete calculation. Column (a) corresponds to the case where the final states 
are pure Hartree-Fock ones. These final states can be renormalized by means of energy and 
spectroscopic factors, either calculated with the same Skyrme force, and we have these values 
in the SIII case in p2| so we can show the corresponding results for the branching ratios in 
column (b), or taken from empirical estimates ||26|| , and the branching ratios obtained in this 
way are shown in column (c). These calculated results are compared with the experimental 
values quoted in Ref. [|]]. The agreement of the total sum of branching ratios with the 
experimental value in the SIII case and in column (b) is quite remarkable, as this is the 
most consistent fashion to perform a theoretical calculation. Minor discrepancies are found 
by looking in the same column at individual channels. On the other hand, calculated decay 
branching ratios depend sensitively on the Skyrme force employed, even when predicted 
peak energy and strength are practically the same. By examining the IAR wave function 
obtained from our calculation, we have noticed that the ph amplitudes corresponding to a 
change in the principal quantum number n are not negligible and this indicates that a simple 
picture of the IAR as given by a T_|0) wave function is slightly too schematic, at least in 
the nucleus we have considered. 



V. GAMOW-TELLER RESONANCE RESULTS 

We follow a similar pattern in presenting the results for the GTR in 208 Bi. Figs. ^ and 
|5] show the strength functions calculated respectively with SIII and SGII interactions and 
Table |I| includes centroids, widths and percentage of strength of these distributions. With 
SIII, at least two main states contribute: besides the main bump at about 21.50 MeV, a 
smaller one appears at lower energy. In the case of SGII, there are two smaller bumps 
in addition to the main one at about 23 MeV. One can understand this broadened line 
shape by looking at the underlying complex states coming out of the diagonalization of the 
effective Hamiltonian Q2.2| ). There are several important complex states spread over the 
energy interval 18-25 MeV. This is in contrast with the IAR case and it can be understood 
by relating to the larger density of T — 1 states compared to that of T states, in this 
energy region of the 208 Bi spectrum. The effect of coupling is much stronger for GTR 
and consequently it admixes TDA states more than for IAR. The percentage of the total 
strength of the GTR operator /3_ = \ Z)u=o ±i ^2t=i cr ^(*) r -(*) (the so-called Ikeda sum rule, 
equal to 3(N — Z)), is about 61% in the energy region 18-24 MeV for SIII, and about 68% 
in the interval 19.5-24.5 MeV, for SGII. This contrasts with the calculation with continuum 
only, where about the same fraction of strength is exhausted by the main narrower peak. In 
the complete calculation the remaining strength is found outside the forementioned interval 
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(E < , E>) and is rather fragmented. Experimentally, the strength concentrated in the GTR 
peak is about 50% of the total but this percentage rises up to about 60-70% if the whole 
energy region with some strength around the main bump is considered P7j . All the present 



theoretical results were calculated with an averaging parameter A equal to 250 keV. Again, 
the stability of the results has been checked by varying A. 

Besides the peak energy, we can extract from the strength distribution a mean energy, 
given by 

fjsZ du uS(u) 

(v) = — I , (5.1) 

!§> du> S(u) 

and this turns out to be 21.11 MeV (SIII) or 22.43 MeV (SGII). Both values of (u) are 
somewhat larger than the experimental excitation energy of GTR which is 15.6 MeV in 
208 Bi, corresponding to 19.2 MeV with respect to the 208 Pb ground state!. As for the width, 
the FWHM is not well defined in this case due to the double or triple peak. If we extract a 
value for the variance a from the strength distribution, by defining 

^ 2 = fj^du (u-(u>)) 2 S(u>) 
Xb< duJ S ( u ) 

we obtain a full width in Gaussian approximation (T =2. 4a) of about 3 MeV for both 
interactions (slightly larger for SGII), to be compared with an experimental value of 3.8 MeV 
@]. Of course, a model in which only the simplest class of "doorway states" is considered 
and not their full hierarchy, cannot but underestimate the width. The same model seems to 
predict correctly the giant monopole width || which is mainly due to Landau damping. 

The overestimation of the peak energy, and the missing part of the damping effect, is 
probably the cause of a result for the branching ratios which correspond to a higher yield 
than the experimental one, as shown in Tables [V] and [VT[ Here, as the strength is spread over 
several MeV, we cannot calculate the branching ratios at a definite energy, as it was done 
for IAR, and we rather make an average over the energy interval (E < , E>) of the numerator 
and denominator of Q2.ll ): 



S C (GTR) = 7 ^ Y . (5.3) 

Tables [V] and |V| are displayed similarly to Tables fTJ and [V] except for the column (b) of 
the theoretical findings which we discuss below. First, one can notice that the discrepancy 
with experiment is reduced when going from self-consistent continuum RPA (column "only 
W^") to the complete model. The results of SGII are in reasonable agreement with the data. 
We must also mention that preliminary indications from the MSU experiment J7J leave open 
the possibility that experimental branching ratios might be larger than found in fl. On the 



3 In order to reach this value, we have to add to the 15.6 MeV the mass difference between the two 
Bi and Pb isotopes, as well as the mass difference between neutron and proton which is missing 
when the GTR is built on the basis of ph excitation. 
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other hand, the sensitivity of predictions with respect to the choice of the interaction is still 
present. We also show in column (b) of Tables and [Vj what the proton yield would be, 
if the excitation energy corresponded to the experimental one (19.2 MeV). The branching 
ratios of column (b) have indeed been obtained by calculating the escape amplitudes (|2.13|) 
with the experimental outgoing proton energies e. The empirical spectroscopic factors pEfl 
already used in the case of IAR are also included. This renormalization obviously decreases 
the branching ratios and they are in good agreement with experiment in the case of SGII 
force. The final difference between theory and experiment, in the case of SIII force, is of the 
same order, of even smaller than the one found in the case of the giant monopole resonance 
||, and simply indicates once more the limit one could expect in this kind of detailed study, 
due to the present uncertainity in the nuclear dynamics at low energy. 

We can finally compare our results obtained without the spreading terms with those 
of Muraviev and Urin |J which are calculated in a continuum-RPA approach, where phe- 
nomenological inputs for the mean field and residual interaction are used. The large differ- 
ences between the two sets of results cannot be explained by invoking the approximation 
made in the present work in treating the single-particle continuum, as already discussed after 
Eq. ( |2.3| ). It should be attributed to the sensitivity of branching ratios to the details of the 
inputs. As we have seen here differences in results coming from the same model but using 
two different Skyrme interactions, it is not too surprising that results from phenomenolog- 
ical RPA calculations, based on a Woods-Saxon potential and a residual interaction of the 
Landau- Migdal type can also differ. A similar situation also exists in the giant monopole 
resonance case |28| . 



VI. CONCLUSIONS 

Within linear response theory, a microscopic model of collective excitations based on 
RPA states coupled to doorway states composed of lp-lh configurations plus a collective 
vibration, and to lp-lh continuum states has been applied to the study of the properties 
of two important charge-exchange modes, the IAR and GTR in 208 Bi. The two types of 
coupling are intended to describe the essential physical mechanisms leading to the spreading 
of the collective mode and its decay by nucleon emission. The model contains no free 
parameter but depends on the choice of the effective nucleon-nucleon interaction. 

The physical quantities which can be calculated are the strength distributions and the 
cross sections of nucleon decay into various channels. From these quantities centroid energies, 
widths and particle decay branching ratios are obtained. The calculated results show some 
sensitivity to the choice of the effective interaction. 

In the IAR case, both SIII and SGII interactions predict the correct peak energy, and 
SIII gives a good description of the data for the decay branching ratios whereas SGII tends 
to underestimate the branching ratios as well as the total width. 

For the GTR, interactions SIII and SGII overestimate the mean energy respectively by 
2 and 3 MeV, while they both predict a total width of about 3 MeV, i.e., 75% of the 
experimental value. Looking at branching ratios, it appears that SGII comes very close to 
experiment if one corrects for the overestimate of the GTR energy. The discrepancy with 
experimental values (in the case of SIII) is of the same order or even smaller than that 
already found in the case of the isoscalar monopole resonance using the same model. This 
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level of agreement with experiment indicates the present limits of existing Skyrme effective 
interactions. 
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APPENDIX A: MATRIX ELEMENTS OF THE EFFECTIVE HAMILTONIAN 

We show in this Appendix how to evaluate matrix elements of the escape and spreading 
terms and of the effective Hamiltonian (|2.2|), on a RPA or TDA basis. We consider 



explicitly the RPA CctSGj clS the reduction to TDA comes out quite transparently. First we 
suppose that the effective Hamiltonian contains only the RPA and the escape term, and we 
derive the explicit form of Eq. (|2.5|) in this case. Then, we show how the spreading term 
can be taken into account. 

We denote a generic element of the RPA basis by \n), its energy being uj n . Another state 
with energy — u n is present in the basis (see chapter 14 of [^|), and we denote it by \n). The 
creation operators of these states are respectively and O^. The creators 0\ of the states 
resulting from the diagonalization of the effective Hamiltonian ( p.2|) , can be expressed 



as a linear combination of the RPA creators, since these are a complete set in Qi, 

Ot = E F n ] Ol ~ F^Ol (Al) 



w„>0 

where the amplitudes F and F are in general complex numbers. 
The eigenvalue equation for the effective Hamiltonian is 

n,oi]^(n v -i^-)ot, (A2) 

where we have fixed at a definite value the energy at which the effective Hamiltonian must 
be evaluated, and omitted to indicate explicitly this energy dependence. Using the property 

[Q.HQ^O^^^Ol (A3) 
(a minus sign appears in the r.h.s. if we consider OjJ, we obtain 

E F^\u n Ol + [W\ OH) - F^(-u n Ol + [W\ Ol\) = (n„ - £ F^O{ - F^O{. 

o; n >0 u> n >0 

(A4) 

If we take the expectation values of (|A4[), first between (0\O m and |0), then between (0\O m 
and |0), we build a set of 2n> equations (n> is the number of positive RPA eigenvalues) 
whose matrix form is the one of (|2T 
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In ( |2.5| ) each submatrix has dimension n > . The submatrix T> is diagonal and contains 
the contribution of the term Q 1 HQi of ( |2.2|) , 



The contributions to the submatrices Ai coming from the escape term are given by 

(A 1 )Z c n = (0\[O m ,[W\Ol]\ |0>, 
(A 2 )Z c n = -(o\ [0 m ,[W\0l]\ |0>, 
(As)~ = (o\[d m ,\w\oH\ |0), 

(A,)Zn = -(0\\O m ,[W\Ot]\\0). (A6) 



To work out explicitly these matrix elements we approximate W' with expression ( [2.3|) (see 
the discussion in Sec. [TI] about this point). The operator K which appears in Q2.3| ) satisfies 
(as said in Sec. |T|but see also [fL3[]) the following equation, 

K 'Qi TT~ ~Qi = Qi- (A7) 

The propagator UJ _^ 0+ie can be expanded in its partial wave components labelled by 
and each component gij(ri,r2) is written by means of the regular and irregular solutions 
(fij(r) and gij(r)) of the corresponding radial Schrodinger equation with energy u, 

2m* 

9ij{ri, r 2 ) = -— 2-/y(r<)0y(r>) • W~\ (A8) 
n. 

where m* is the nucleon effective mass, r> (r<) are the larger (the smaller) between n and 
r2, and is the Wronskian of /y and gij [[SO]- In this form the propagator is more easily 
inverted to obtain K. This is a ph operator diagonal in the hole and for a given hole h we 
label as Kh the corresponding particle operator. 

The matrix elements of the escape term are then written in terms of the ones of Kh and 
of the RPA amplitudes xj$ and Y^' as 

(AOmn = ~^2( £ p ~ £ h){xfflX$ + Y^Y^) - ^2 (p'\ K h\p}Shh'(Xp7l l }x^ l ) + Y$$Yffi), 

ph ph,p'h' 

C4,)~ = E [(% " Zh)5pp>6hh> + (jp^K h \p)5 hh ]{X$Y$ + Y$X$), 

ph,p'h' 
(^Imn (• / ^2) mrl ; 

= (-^Omrr (A9) 

We now turn to the calculation of the spreading term. It is known (see, e.g., the first chapter 
of [|TTJ) that diagonalizing an effective Hamiltonian including this term besides the RPA one, 
and with a set of 2p-2h configurations in the Q 2 space, amounts to solve 
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where A and B are the RPA matrices and has the matrix elements defined by fl2.4|) . 



Eq. ( A10 ) is the eigenvalue problem in the SRPA scheme mentioned in the text. A similar 
pattern holds even if a different choice of Q 2 is madei, like in the present case, and therefore, 
in order to obtain the spreading matrix elements which must be inserted in ( |2.5|) one has to 
transform the spreading part of Eq. ( |A10| ) which is written on the ph basis of Qi, and write 



it on the RPA basis. The result is the spreading contribution to ( p.5|) , which reads 
(A)mn = E ^ph^'h'i^^ph X^li + W^ p , h ,(—uj)Y p ™ Yj^,, 

ph,p'h' 

( A ^spr _ w i / A y{m) v {n) j,,;* / \y(m) y{n) 

\^2)mn — vv phtfh'V^ )^ p h I p'h' ' vv ph,p'h'\ u ) 1 ph ^p'h'i 

ph,p'h' 

(^)mn — E ^ph,p'h'{° J ) Y P h' Xp? h , +W^ p , h ,(-u)X p ™ Y^, 
ph,p'h' 

(A^) S mn = E Wp hip t h ,(u))Y ph Y p , h , + Wp hplh ,(— Uj)Xp h ^X p ,y. (All) 
ph,p'h' 

The matrix elements W php , h , can be evaluated as discussed in Sec. fTl]. We give here their 
final expressions, 



k=l 



L,n,p" W ~ K + £ P" ~ £ h) + iA fp 

dr x u v > {r x )u p » {r^vir^g^ (r x ) / dr 3 u p (r 3 )u p » (r 3 )v(r 3 )g^ (r 3 ) , 



^(2) = s P:p/ s jhJh/ ^ 



f_l)i+V+^'+ L+J f ,\ / 1 

1^(3) = E 7 ^ f u , f f r (jp||niiv)^iiybibfc> 

^ uj - (u n + e p - £ h t) + 2A [ Jh' J P > L> J 
• J dri«p(ri)M p /(ri)t;(ri)^ L) (ri) ^ rfr 4 M h (r 4 )'u^(r4)i)(r4)^ i) (r4), 

^(4) = E 7-T ■ : T^-X ? • r P II^HV)(^II^||J,) 



L,n 



{uj n + £p/ - e h ) +iA { jh' j P ' L 



■ J dr 2 w p (r 2 )iy(r 2 )v(r 2 )^ L) (r 2 ) J dr 3 u h (r 3 )u h ,(r 3 )v(r 3 )g^\r 3 ), (A12) 
where the symbol j denotes 2 j + 1. 



4 The eigenvalue problem has the same form but not all the statements about SRPA contained 
e.g. in p3 are, strictly speaking, still valid. 
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APPENDIX B: DECAY BRANCHING RATIOS 



In (T3| the cross sections a c and a exc (see Sec. [TI|) are defined as 



w 

<W = ~l^rIrnJ2 T^F* F T ) UU ,(^ ~ 0/ - ^f)' 1 ^, (Bl) 

where all the quantities are defined in the text, apart from k which is here the wave vector 
corresponding to the initial state of the projectile and the elements of the T-matrix, in which 
one state is the initial state of target (labeled by 0). These matrix elements are (see also 

0) 

T v0 ={ X i f ) ^\V\% X ^ ) ), (B2) 

where the initial and final wave functions of the projectile and ejectile are indicated by x 
and those of the target nucleus by Here, V is an effective projectile-target interaction. 
We consider for our present purposes its r • r (a ■ or ■ f) channel for the case of IAR (GTR) 
excitation. Either of these operators is labelled by O, so the interaction is Voif\ . . .ta\ r p ) ■ O 
where f p is the projectile spatial coordinate. We are in the present work not interested in 
the reaction mechanism, provided it can excite the target in a state with definite quantum 
numbers. Therefore, we can make simplifications about this mechanism and we use standard 
PWBA methods to evaluate the element (|B~^). Plane wave functions are substituted to the 
initial and final wave functions of the projectile and ejectile, and an effective factor K takes 
care of the reduction of the reaction probability due to distortion which would come from a 
realistic estimate made by using an optical potential. In a straightforward way (|B2|) becomes 



T v0 = KV o (q)(u\e^-O\0), (B3) 

where <f is the momentum transferred by the projectile and V(q) is the Fourier transform of 
the interaction V. The exponential function can be decomposed as usual, 



v^EHjV^A + l- 3x (q-r)Y X0 (r), (B4) 



e 

A=0 



using Bessel functions and spherical harmonics. If we restrict to the A-value under consid- 
eration (0 as both IAR and GTR have no change of angular momentum with respect to 
the ground state in the space part of their wave function) and truncate the expansion of 
the Bessel function at the lowest order (low momentum transfer approximation), we obtain 



from (|B4| ) a constant which can be included in the definition of V. Finally we are left with 



the approximate expression 

T u0 ~ constant ■ (u\O\0). (B5) 

If we use this expression it is easy to show that the branching ratio, i.e. the ratio between 
the two cross sections defined in (|B1| ), is given by (|2.11|) where 
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S vv > = (/v|O|0)(z/ / |O|0)*, (B6) 

and the other quantities are defined in Sec. ||. 

We conclude this Appendix by proving that the branching ratios defined in ( [2.1 1| ) are 
real. We show that the complex conjugate, 

r* _ -■y itu' ,cStu'[i. UJ y' ~ Uv) ~ I (]V + r ^)] I sr^i t?* r?V c* ^ / R7 \ 

is equal to S c . Using the fact that the matrix F is complex orthogonal one obtains 

(F*F T y uul = {F*F T ) V , U . (B8) 



We can insert in ( p7|) this last expression, together with the analogous properties of S and 
7, S* u , = S v i v and 7*^/ c = 7 I /„, C . Exchanging the indices v and 1/ completes the proof. 
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FIGURES 



FIG. 1. Diagrammatic representation of the four terms whose sum gives the matrix element 
. The analytic expressions are shown in Appendix A, Eq. ( |A12 ). 

FIG. 2. Strength distribution of IAR in 208 Bi calculated with the interaction SIII. The full line 
refers to the complete calculation, whereas the dashed line gives the result with only the continuum 
coupling. 

FIG. 3. Strength distribution of IAR in 208 Bi calculated with the interaction SGII and displayed 
as in Fig. ||. 

FIG. 4. Strength distribution of GTR in 208 Bi calculated with the interaction SIII. The two 
lines have the same meaning as in the case of IAR. 

FIG. 5. Strength distribution of GTR in 208 Bi calculated with the interaction SGII and dis- 
played as in Fig. f|. 
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TABLES 



TABLE I. Vibrational states of Pb, resulting from a coordinate-space RPA calculation em- 
ploying the SGII force, with energy below ~ 20 MeV, and exhausted percentage of EWSR larger 
than 1%. They are used to build the "doorway states", that is, the configurations which are more 
complicated than the lp-lh configurations and therefore play a role in the damping process of giant 
resonances. The experimental values in this Table are taken from |23| . 

Theory Experiment 



3 n Energy B(EA) Percentage Energy B(EA) 

[MeV] [spu] of EWSR [MeV] [spu] 

(or perc. of EWSR) 



2+ 


4.78 


3.28 


7 


4.085 


6.2 




5.13 


4.36 


10 


4.923 












5.036 












5.128 






1 1 i 




l o 


1 n r 


1 /OFj VV OIL — / U 1 


3~ 


2.64 


37.68 


25 


2.614 


34.2 




13.6 


0.62 


1.6 








14 


0.37 


1 








16.53 


0.41 


1.3 








17.55 


0.30 


1 








18.17 


0.40 


1.4 








20.68 


10.13 


40 


20.9 


(%EWSR=36) 


4+ 


4.82 


2.52 


1.1 


4.323 


23 




5.71 


11.21 


5.8 


5.690 






12.05 


9.61 


10.5 


12. 


(%EWSR=10) 




12.44 


4.52 


5.1 








13.17 


2.43 


2.9 








14.66 


1.43 


1.9 








15.14 


0.80 


1.1 






5" 


3.58 


16.37 


2.6 


3.198 


17.4 
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TABLE II. Averaged quantities extracted from the strength distributions of IAR and GTR. In 
the case of IAR only one peak appears, as expected, and we display the energy, the full width and 
the percentage of strength under the peak. In the case of GTR we take into account the whole 
energy region (E K , E>) in which bumps of the strength function are visible. The values of E < and 
E >} as well as the definitions of the mean energy and width, are given in the text (Sec. 0). 



Isobaric analog resonance 








Theory 










SIII 




SGII 


Experiment |25|] 


Mean Energy 




18.49 MeV 




18.64 MeV 


18.8 MeV 


Width 




152 keV 




99 keV 


232 keV 


Percentage of strenj 




97% 




97% 


~ 100% 


Gamow- Teller resonance 








Theory 










SIII 




SGII 


Experiment Q 


Mean Energy 




21.11 MeV 




22.43 MeV 


19.2 MeV 


Width 




3 MeV 




3.1 MeV 


3.7 MeV 


Percentage of strenj 


$th 


61% 




68% 


~ 60-70% |27| 



TABLE III. Branching ratios for the proton decay of IAR in 08 Bi to neutron valence hole states 
of 207 Pb, obtained by using the SIII force. In column (a) the final state is a pure HF configuration, 
whereas in other columns its energy and wave function are corrected either by means of a consistent 
SIII calculation (b) or by using empirical values (c). 



Theory 



Decay 


only W T 




W T + w l 




Experiment 


channel 




(a) 


(b) 


(c) 


1 


Pi 


0.472 


0.346 


0.253 


0.237 


0.22±0.02 


2 

P3 

2 


0.396 


0.287 


0.238 


0.196 


0.34±0.04 


i i3 


0.015 


0.011 


0.008 


0.010 




2 
f 5 


0.117 


0.086 


0.065 


0.061 


includ. in p3 


2 

fr 


< 10~ 3 


< 10-3 


< 10-3 


< lO- 3 


0.015±0.007 


2 

hs 


< 10~ 3 


< lO- 3 


< 10-3 


< 10-3 




EcBc 


1. 


0.730 


0.564 


0.504 


0.575±0.07 
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TABLE IV. Same as Table 



III. 



, in the case of the SGII force. 



Theory 



Decay 


only W 1 


TT rT 1 ~\T T 1 

W 1 + W l 




Experiment 


channel 




(a) 


(c) 




D i 

2 


0.448 


0.171 


0.137 


0.22±0.02 


P | 


0.514 


0.198 


0.157 


0.34±0.04 


i 13 

2 


0.016 


0.007 


0.006 




fs 

2 


0.016 


0.008 


- 0.006 


includ. in p3 


fr 

2 

ll9 


0.006 


0.002 


0.004 


0.015±0.007 


< 10~ 3 


< 10~ 3 


< 10-3 






1. 


0.386 


0.310 


0.575±0.07 



TABLE V. Branching ratios for the proton decay of GTR in 208 Bi, obtained by using the SIII 
force. In column (a) the final state is a HF configuration. In column (b) we present the results of a 
calculation where the energy of the final neutron hole state is corrected in such a way the escaping 
particle has the same energy as in the experiment Q, and spectroscopic factors [26| are used. 







Theory 






Decay 


only W T 


W T + wJ- 




Experiment 


channel 




(a) 


(b) 


1 


Pi 


0.265 


0.037 


0.022 


0.013±0.002 


2 

PI 


0.432 


0.055 


0.033 


0.023±0.003 


2 

ii3 

2 


0.009 


0.001 


0.001 


0.002±0.002 


fs 

2 

fr 

2 


0.278 


0.051 


0.030 


includ. in p3 


0.011 


0.009 


0.005 


0.003±0.002 


h g 


0.005 


0.001 


0.001 




Ec^c 


1. 


0.154 


0.092 


0.041±0.009 
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TABLE VI. Same as Table [V], in the case of SGII force. 



Theory 



Decay 


only W 1 




W 1 + W + 


Experiment 


channel 




(a) 


(b) 


IL_|I 


Pi 

2 


0.223 


0.033 


0.018 


0.013±0.002 


pa 


0.418 


0.035 


0.019 


0.023±0.003 


2 


0.014 


0.003 


0.001 


0.002±0.002 


u 

2 


0.319 


0.013 


0.007 


includ. in p3 


fr 

2 

hg 


0.016 


0.010 


0.003 


0.003±0.002 


0.010 


0.001 


< 10~ 3 




EcBc 


1. 


0.095 


0.048 


0.041±0.009 
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